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We study the structure of hadronic protoneutron stars within the finite tempera- 
ture Brueckner-Bethe-Goldstone theoretical approach. Assuming beta-equilibrated 
nuclear matter with nucleons and leptons in the stellar core, with isothermal or 
isentropic profile, we show that particle populations and equation of state are very 
similar. As far as the maximum mass is concerned, we find that its value turns out 
to be almost independent on T, while a slight decrease is observed in the isentropic 
case, due to the enhanced proton fraction in the high density range. 



I. INTRODUCTION 



A protoneutron star (PNS) is formed after a successful supernova explosion and consti- 
tutes for several tens of seconds a transitional state to either a neutron star or a black hole 
l|. Initially, the PNS is optically thick to neutrinos, that is, they are temporarily trapped 
within the star. The subsequent evolution of the PNS is dominated by neutrino diffusion, 
which first results in deleptonization and subsequently in cooling. After a much longer time, 
photon emission competes with neutrino emission in neutron star cooling. 

In this paper, we will focus upon the essential ingredient that governs the macrophys- 
ical evolution of neutron stars, i.e., the equation of state (EOS) of dense matter at finite 
temperature. We have developed a microscopic EOS in the framework of the Brueckner- 
Bethe-Goldstone (BBG) many-body approach extended to finite temperature. This EOS 
has been successfully applied to the study of the limiting temperature in nuclei 2l, IS]. The 
scope of this work is to present results on the composition and structure of these newly born 
stars with the EOS previously mentioned. 

Two new effects have to be considered for a PNS. First, thermal effects which result in 
entropy production, with values of a few units per baryon, and temperatures up to 30-40 
MeV jj, sj]. Second, the fact that neutrinos are trapped in the star, which means that the 
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neutrino chemical potential is non-zero. This alters the chemical equilibrium and leads to 
compositional changes. Both effects may result in observable consequences in the neutrino 
signature from a supernova and may also play an important role in determining whether or 
not a given supernova ultimately produces a cold neutron star or a black hole. 

This paper is organized as follows. In Section 2 we briefly illustrate the BBG many-body 
theory at finite temperature. Section 3 is devoted to the study of the stellar matter composi- 
tion at constant temperature or entropy, and the resulting EOS. In Section 4 we discuss our 
results regarding the structure of (proto) neutron stars, in particular their maximum mass. 
Finally, conclusions are drawn in Section 5. 



II. THE BBG THEORY AT FINITE TEMPERATURE 



In the recent years, the Brueckner-Bethe-Goldstone perturbative theory 



extended in a fully microscopic way to the finite temperature case [2], according to the 
formalism developed by Bloch and De Dominicis In this approach the essential ingredient 
is the two-body scattering matrix K, which, along with the single-particle potential U, 
satisfies the self-consistent equations 

+ Rej2(hh\v\k>k' A; n ^m-n(K)} {mKmhh) (1) 

and 

U{h) = J2 n (k 2 )(k 1 k 2 \K(W)\k 1 k 2 ) A , (2) 

fc 2 

where ki generally denote momentum, spin, and isospin. Here V is the two-body interaction, 
W = + Ek 2 represents the starting energy, E\~ = k 2 /2m + U(k) the single-particle 
energy, and n(k) is the finite temperature Fermi distribution. We chose the Argonne Vi§ 
nucleon-nucleon potential [8j as two-body interaction V, supplemented by three-body forces 



6| has been 



(TBF) among nucleons, in order to reproduce correctly the nuclear matter saturation point 
On ~ 0.17 fm~ 3 , E/A ~ —16 MeV. We have adopted the phenomenological Urbana model 
[91], which consists of an attractive term due to two-pion exchange with excitation of an 
intermediate A resonance, and a repulsive phenomenological central term. In the BBG 



approach, the TBF is reducec 
position of the third particle 



to a density- dependent two-body force by averaging over the 
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In order to simplify the numerical procedure, we introduce the so-called Frozen Correla- 
tions Approximation, i.e., the correlations at T 7^ are assumed to be essentially the same 
as at T = 0. This means that the single-particle potential Ui(k) for the component i can 
be approximated by the one calculated at T = 0. The accuracy of that approximation has 



been checked in Ref. {2]. 

Within this approximation, for a fixed density p = n(k) and temperature T, we solve 
self-consistently Eqs. ([1]) and along with the equation for the free energy density, which 
has the following simplified expression 



i 

where 



(3) 



s l = - ( n i( k ) ^ni(k) + [1 - n i{k)} ln[l - m(k)]) (4) 
k 

is the entropy density for component % treated as a free gas with spectrum Ei(k). Finally, 
the pressure can be easily calculated as the derivative of free energy with respect to the 
density, keeping fixed the temperature, i.e., 

III. COMPOSITION AND EOS OF HOT STELLAR MATTER 

For hot and neutrino trapped stellar matter containing only nucleons as relevant degrees 
of freedom, the composition is determined by the requirements of charge neutrality and beta 
equilibrium, which read explicitly 

^ qiXi + ^2 1 lXl = ' ( 6 ) 

i I 

fM n -fM p = p; e - /v = /v + . (7) 

In the expression above, xi = Pi/pB represents the baryon fraction for the species i, pB the 
baryon density, and qi the electric charge. The same holds true for the leptons, labelled by 
the subscript I. It turns out that the dependence on the proton fraction can be approximated 



by the so-called parabolic approximation which reads for the free energy as 

j(p,x p ,T) « j(p,x p = 0.5, T) + (1 - 2x p ) 2 F sym (p,T) (8) 
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being F sym the symmetry energy 



F F 

F sym (p,T) « —(p,x p = 0,T) - — (p, x p = 0.5, T) 



(9) 



In the parabolic approximation, the nucleon chemical potentials can be easily calculated as 



/i n (p, x p , T) 



df_ 

dp n 



1+P 



fi p (p,x p: T) = ^- 
op p 
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(10) 
(11) 



whereas the chemical potentials of the noninteracting leptons are obtained by solving numer 



ically the free Fermi gas model at finite temperature. Further details are given in Ref. 
Because of trapping, the numbers of leptons per baryon of each flavor I = e, p, 



Yi — Xi — Xi + Xj, 



X 



(12) 



are conserved on dynamical time scales. Gravitational collapse calculations of the white- 
dwarf core of massive stars indicate that at the onset of trapping the electron lepton number 
is Y e w 0.4, the precise value depending on the efficiency of electron capture reactions during 
the initial collapse stage. Moreover, since no muons are present when neutrinos become 
trapped, the constraint = can be imposed. We fix the Yi at those values in the 
calculations for neutrino-trapped matter. 



A. Isentropic description of hot stellar matter 

As known from dynamical simulations, thermal effects result in an approximately uniform 
entropy per baryon of S/A ~ 1—2 across the star [4]. Therefore, within the same theoretical 
framework, it may be useful to switch from an isothermal to an isentropic description of hot 
and trapped stellar matter, and compare the respective outcomes. In our approach, we 
proceed as follows. Once we determine the composition and the free energy per particle 
F I A at a given density p, for several values of the temperature T, we calculate the entropy 
per baryon (in units of the Boltzmann constant) from the thermodynamical relationship 

Hence, for fixed S, we find how the temperature changes as a function of the density p, as 
shown in Fig. [TJ The lower curve displays the calculation for S/A — 1, whereas the upper 
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curve is the calculation for S/A = 2. We notice that the temperature is a monotonically 
increasing function of the density, and that the reachable values are larger for larger values 
of the entropy. We remark that the range of temperatures coincides with the values chosen 
for the isothermal description of a protoneutron star 12 1. 

Following this procedure, we obtain the stellar composition and the corresponding free 
energy at a given nucleon density p. Finally, we can calculate the total energy E/A from 
the thermodynamical relationship 



E F S 

— (p, x p , S) = — (p, x p , S) + T(p, Xp)- , (14) 



and the equation of state is 



The pressure vs. density relationship is the fundamental input for calculating stable con- 
figurations of a protoneutron star. These can be obtained from the well-known hydrostatic 
equilibrium equations of Tolman, Oppenheimer, and Volkov [13( for the pressure P and the 
enclosed mass m 

dP(r) = Gm(r)e(r) [1 + + 

^ r r 2 , _ 2Gm(r) ' V 1 "/ 

r 

d ^P- = We(r) , (17) 
dr 

once the EOS P(e) is specified, being e the total energy density (G is the gravitational 
constant). For a chosen central value of the energy density, the numerical integration of 
Eqs. ffTBT) and ffT7|) provides the mass-radius relation. For simplicity, we have attached a 
cold crust to the stellar core, by joining the hadronic equations of state above described 
with the ones by Negele and Vautherin [ijj] in the medium- densit y r egime (0.001 fm~ 3 < 



p < 0.08 fm~ 3 ), and the ones by Feynman, Metropolis, and Teller [15( and Baym, Pethick, 
and Sutherland [16| for the outer crust (p < 0.001 fm~ 3 ). 

A more realistic treatment of the transition from the hot interior to the cold outer part 
has a dramatic influence on the mass-central density relation in the region of low central 
density and low stellar masses, as discussed in Ref . 17| . However, the maximum mass region 
is not affected by the structure of this low-density transition region. 
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IV. RESULTS AND DISCUSSION 

Let us now discuss first the populations of beta-equilibrated stellar matter, by solving 
the chemical equilibrium conditions given by Eqs. ((Zj), supplemented by electrical charge 
neutrality, and baryon and lepton number conservation. 

In Fig. [2] we show the particle fractions as a function of the nucleon density, for different 
values of the temperature, respectively T=20,40 MeV (upper panels), and entropy S/A— 1,2 
(lower panels). We do not observe sizeable differences between the two calculations, except 
for the neutrino fraction at extremely low densities, and the muon percentage in the high 
temperature/entropy case. 

The equation of state for beta-stable and neutrino trapped matter is displayed in Fig. El 
where the pressure is reported as a function of the nucleon density. We observe that the 
equation of state stiffens with increasing T or S mostly in the medium-low density range. 
On the contrary, at large density thermal effects play a minor role, and curves of different 
T or S show a quite similar behavior, as already found in Ref. [l^ . 

Finally, in Fig. H] we display the gravitational mass as function of the radius (left panels), 
and the central energy density e normalized with respect to the saturation value eo (right 
panels). The dashed (dot-dashed) line represents the calculations for the low (high) temper- 
ature or entropy case. We find that the values of the maximum mass are quite stable and 
almost independent on S or T. In particular, we observe that the value of the maximum 
mass slightly decreases with increasing S. The reason is that in hot and trapped stellar 
matter, the proton fraction increases with increasing S or T at large values of the nucleon 
density, and this leads to a softening of the equation of state, with consequent decrease of 
the maximum mass. This result is in contrast with the findings of Prakash et al. [1( 5 where 
the critical mass increases with increasing entropy. This can be due to the fact that in 



Ref. |l| the interaction is mostly local, with only a non-local correction, and therefore only 
the kinetic part contains a dependence on temperature. In our calculations the whole inter- 
action part is temperature-dependent, due to the intrinsic non- locality of the single-particle 
potentials, and therefore the temperature dependence can be quite different. 

In the same figure we also show the results for cold neutrino-free neutron stars (thin solid 
line). One notes in particular that the maximum masses of the protostars are always slightly 
smaller than that of the neutron star, which excludes a delayed collapse while cooling down 
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V. CONCLUSIONS 

We have modeled nucleonic protoneutron stars using a realistic BBG equation of state and 
assuming either isothermal or isentropic conditions. Heavy protostars are fairly insensitive 
to that choice and their maximum masses are slightly smaller than that of the cold neutron 
star. 
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FIG. 1: Temperature vs. nucleoli density for neutrino-trapped matter with entropy S/A=l (dashed 
line) and S/A=2 (solid line). 
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FIG. 2: Relative populations for neutrino-trapped matter as a function of the nucleon density. 
Calculations are shown for the isothermal case (upper panels), and for the isentropic one (lower 




FIG. 3: Pressure as a function of the nucleon density for the isothermal case (left panel) and the 
isentropic one (right panel). 




FIG. 4: Mass-radius (left panels) and mass-central density (right panels) relations for isothermal 
(upper panels) and isentropic (lower panels) conditions. 



